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Abstract 

We discuss a dispersion relation in the photon mass and show how (in prin- 
ciple) model-independent constraints on the parton distribution functions of 
the photon, notably a momentum sumrule, can be obtained. We present two 
sets of parametrizations, SaS 1 and 2, corresponding to two rather extreme 
realizations of the non-perturbative part. Inclusive electron scattering off 
a real photon is found to be insufficient to constrain the non-perturbative 
components. The additional sensitivity provided by the photon virtual- 
ly is outlined. Previous approaches to model the non-perturbative input 
distributions are commented upon. 
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1. Introduction 

Perturbative QCD predicts only the Q 2 evolution of the parton distribution functions 
(PDFs) of the photon fi(x,Q 2 ) via a set of inhomogeneous differential equations of 
the first kind. Hence the solutions fi{x, Q 2 ) require the specification of the PDFs at 
some Q 2 = Ql- Two ways exist to determine these non-perturbative input distributions 
fi(x,Qo)- The first one is analogous to the determination of hadronic PDFs: At Qo 
large enough to be safely within the perturbative regime (Qo ~ 2GeV), the parameters 
of the input distributions f2( x ,Qo) (shapes and normalizations) are fitted |I| to the ex- 
perimentally measured distributions involving the PDFs of the photon, which thus far 
means to the available F^x^Q 2 ) data. Since these data are currently restricted to large 
x, only the n-valence distribution is known with some confidence. In particular, basically 
no constraint on the gluon distribution of the photon exists today in such an approach. 

In the second approach one pretends to know the input distributions at some very 
low scale Qo ~ 0.5 GeV apart from a single, adjustable parameter. The expectation is 
that, at such low scales, the photon should essentially behave like a hadron and, corre- 
spondingly, the PDF of the photon could be identified with appropriate hadronic ones. 
The experimental evidence for this ansatz is, however, rather weak: the only data for Q 2 
below 2 GeV 2 come from the TPC/27 measurement @ of F^ix, Q 2 ) (at an average Q 2 of 
about 0.7 GeV 2 ) and consist of no more than a handful of points in a limited x-range; the 
largest x-bins are moreover plagued by resonance contributions. Nonetheless, basically 
all recent experimental analyses accept the hadron-like parametrization of F^(x, Qq) of 
TPC/27 as the non-perturbative input. The scale Qo is considered as a free parameter 
and fitted to their data, hence disregarding most of the potential of their own F^ (x, Q 2 ) 
data to extract the non-perturbative part. Rather, these analyses merely quantify how 
compatible the more recent data are with the TPC/27 ansatz. The actual, fitted value 
of Qq is, in fact, not a significant number since it is strongly correlated with the size (and 
shape) of the assumed non-perturbative input. 

This correlation is most easily seen by decomposing the PDFs of the (real, i.e. P 2 = 0) 
photon as follows 

f2(x, Q 2 ) - f^(x, Q 2 ) = f:> PT (x, Q 2 , Ql) + f^ p (x, Q 2 , Q 2 ) , (1) 

where the second term on the LHS describes the (properly normalized Z 3 = 1 + 0(a em )) 
probability distribution of a photon to remain a photon 

P a ^(x } Q 2 )=Z 3 5 a ,5(l-x) . (2) 

Being the solution of an inhomogeneous evolution equation, the PDF of the photon can 
always be written as the sum of two terms, as in (||) where the first term on the RHS is 
a particular solution of the inhomogeneous equation with the boundary condition 

f2> PT (x 1 QlQ 2 ) = 0. (3) 

The second term is a general solution of the corresponding homogeneous evolution equa- 
tion and needs a (non-perturbative) input distribution at Q 2 = Q^: 

r a ' NP (x,Ql,Ql) = f~r(x). (4) 

At this point it should also be pointed out that the perturbative Q 2 evolution is not 
treated correctly in the experimental analyses of the photon structure function. The scale 
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Qo is fitted to the respective data on 

F2{x, Q 2 ) = 2xJ2el [ff\x, Q 2 , Ql) + f^ p (x, Q 2 , Q 2 ) } 
g 

ee F^ T {x, Q 2 , Ql) + i^ NP (x, Q 2 , Ql) , (5) 

using the FKP parametrization || of K7 ,PT (;r, Q 2 , Ql). An error in the Q 2 evolution of 
F^x^Q 2 ) arises because the hadronic part is not evolved with Q 2 but kept fixed at the 
input scale, F%' P (x, Q 2 ,Qo) = F2 (x)[TPC/2j]. There is yet another error: the FKP 
parametrization is based on a valence approximation, and hence fails for x < 0.3 [|J. 

The second approach, namely approximating the photonic input distributions at some 
low Q ~ 0.5 GeV by hadronic ones, has also been pursued in theoretical analyses ||. 
Here the input distributions are identified with those of the pion and Qo is fixed by 
theoretical prejudice. In order to have an adjustable parameter, the overall normalization 
of the input distributions is allowed to vary. (This "iT-factor" is actually fitted to high-Q 2 
(Q 2 Ql) data only. In this way one "only" assumes the leading-twist formula to evolve 
perturbatively down to low scales but not to describe all the physics at low scales.) 

It should be stressed that the estimation of the input distributions fa ,NF ( x , Ql, Ql) by 
the ones of the pion using vector-meson dominance (VMD) and the additive quark model 
involves quite severe assumptions. The PDFs of the qq "bound states" of the photon 
need not be the same as those of real vector mesons. Moreover, the PDFs of the short- 
lived p°-meson may well differ in shape from those of the long-lived pion. In addition, 
not only the shapes may differ. Also the relative normalizations of the various PDFs 
can be different, e.g. down-quark to strange-quark distributions, valence to sea to gluon 
distributions. Finally, implicit (but never justified) in this approach is the assumption of 
a momentum sumrule and a constraint on the number of valence quarks. 

2. The dispersion relation in the photon mass 

It will be shown here how a dispersion relation in the photon mass can be used to obtain (in 
principle) model-independent constraints on the PDFs of the photon, namely a momentum 
sumrule, as well as constraints on the valence distribution and the overall normalization. 

The moments of the photonic PDFs (g(n) = Jq 1 dxx n ~ 1 g(x)) can be represented as a 
dispersion integral in the photon mass a 2 (P 2 is the photon virtuality) || 

P a (n, Q 2 , P 2 ) = / — — p a (n, Q 2 , a 2 ) . (6) 
jo eH + P z 

Rather than describing the dispersion integral as the difference between a "point-like" 
part (contribution from the upper limit) and a "hadronic" part (contribution from the 
lower limit), it is more natural to separate short- distance and long-distance parts by a 
scale Qo, since the weight function p a possesses the scaling- violation pattern typical of 
ordinary hadronic PDFs ||. At large values of a 2 , the 7 — > qq transition can be calculated 
perturbatively. At lower values of a 2 one enters the resonance region: non-perturbative 
Regge poles will contribute signalling the appearance of qq bound states. A general ansatz 
for the weight function is therefore (p' = dp/da 2 ) 

Vm(Qo) / 2 \ V m (Q ) / 2 \ 

Pa (n,Q 2 ,a 2 )= £ A v a (n,Q 2 )6 1-^ + £ B^(n, Q 2 ) 5' 1 - ^ 
v=v \ m vJ v=v \ m vJ 

+ 6 (a 2 - Ql) {a a (n, Q 2 , a 2 ) + (3' a (n, Q 2 , a 2 )} . (7) 
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The coefficients , u a and j3 a can be determined as follows. VMD is known to 

well describe photon-hadron interactions over a wide range of energies, from y/s of a few 
GeV up to the HERA energy (200 GeV). Hence, to very good approximation, one may 
neglect A y a and take 

B?(n, Q\ Ql) = (j-) K v (n, Q\ Q 2 ) . (8) 

In order to obtain a a and (3 a one first notices that for Ql <C P 2 <C Q 2 the resonance 
contributions to (0) are suppressed. Moreover, in this limit the PDF of a virtual photon 
(i.e. the LHS of @) can be calculated within perturbative QCD 0. Then one expresses 
these distributions as an integral of "state" distribution functions f2 ,qq {x,Q 2 ,a 2 ): 

f2(x, Q 2 , p 2 ) = f ^r^EK Q\ a 2 ) , (9) 

J P z G Z7T q 

which obey the standard, homogeneous evolution equations with the boundary condition 

I 

f2**(x, a 2 , a 2 ) = K**(x) ^ (V + (1 - x) 2 ) (6 aq + 6^) . (10) 

Equation (^) yields an expression for the sum a a plus (3' a . The decomposition into a a and 
P a is more difficult. Generalized VMD arguments suggest a a <C P a and hence we arrive 
at the final expression for the PDFs of the virtual photon 

P a (n, Q 2 , P 2 ) = fr p (n, Q 2 , P 2 ) + P a ^(n, Q 2 , P 2 ) 
VM f m 2 \ 2 4 



+ r — a ^ fee?) tfV) • (ll) 

•/q* (a 2 + P 2 ) 2 7T V 9 V v ; 

Equation ([TT|) contains three constraints on the non-perturbative distributions. The first 
two follow from the fact that f2' are ordinary mesonic PDFs. Hence they should respect 
the number of valence quarks 

1 = &(n = 1, Q\ Ql) = [ dx Q 2 , Ql) (12) 

and obey the momentum sumrule 

1 = E Pa V (n = 2, Q 2 , Ql) = £ dx x f2> v (x, Q 2 , Ql) . (13) 



a=q,q,g 



The third constraint on the photonic PDFs follows from the observation that the RHS 
of (|ll|) has to be Qo-independent. This fixes the overall normalization and, in turn, gives 
a momentum sumrule for the photonic PDFs: 



1-Z 3 = £ f'dxx f2(x,Q 2 ,P 2 ) 



a=q,q. 

y% Wv + P 2 ) Pv + k (a* + Pi) 2 n \Y q ) 1 J 
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since the perturbative state distributions f^' m (n,Q 2 ,a 2 ) obey a relation analogous to 



This last equation is, in fact, nothing but the observation that the vacuum fluctuations 
of the photons probed in ey interactions are precisely the ones seen in e + e~ — > hadrons. 
The probability per unit cr 2 of their occurrence is given by the cross section of the latter 
reaction 

W-^V -*w> . (1B) 

\a 2 + P 2 J 47r 2 a em v ' 

Using the narrow-width approximation for the resonance (low-mass) contribution and the 
parton-model result at high masses 

o"tot(e + e~ -> X(ct)) = 
X> 2 a em [£j 5 (a 2 - m 2 v ) + (a 2 - Ql) N c f^e 2 ^ ^= (16) 

the result flT4]) is recovered. 

Once Qo has been determined for a given number of vector mesons, flUD tells us how 
Qo has to be changed in order to compensate the inclusion (or omission) of a vector 
meson. For a given number of included vector mesons, the value of Qo can be determined 
from the continuity requirement of the e + e~ annihilation cross section or, say, the total 
7P cross section ||. 

What remains undetermined in the approach sketched above is the shape of the non- 
perturbative input distributions. In line with the argument that hard processes probe 
short time scales, the contributions from the various vector mesons should be added 
coherently. Also, to rather good approximation, an SL^-symmetric sea distribution s(x) 
can certainly be assumed. Then the non-perturbative input distributions f^ ,NP (x, Ql, Ql) 
are given in terms of three distributions, a valence distribution v(x), a gluon distribution 
g(x), and the sea distribution s(x). These distributions can be determined through a fit 
to the available (real photon) F^x^Q 2 ) data, subjected to the constraints (|Tl~D, (|T2]) and 



(0)- 

Two extreme scenarios can be considered. In the first, VMD is restricted to the 
well-established p°, uo, states. The scale Qo is then known to be Qo m 0.6 GeV from 
an analysis of the 7p total cross section ||. This "low-Qo" fit (SaS set 1 distribution 
functions; for details of the fit see |4]]) is essentially a three-parameter fit, two for the 
shape of the valence distribution and one for the normalization of the sea. The shapes of 
both the gluon and sea distributions are hardly constrained by current K7 data and take 
on values of an educated guess. The main theoretical error of this scenario arises from 
the use of perturbation theory down to rather low values of Q 2 . 

The spirit of the second analysis is opposite: take Qo well within the perturbative 
domain (Q = 2 GeV) at the expense of parametrizing the effects of additional vector 
mesons (besides p°, u, and 4>) by a simple factor K to be fitted to the data 

V=V JV V=p°,u),(f> Jv 

This "high-Qo" fit (SaS set 2 distribution functions [|J) contains two additional parameters 
compared to the low-Qo fit) one parameter characterizing the necessary hard component 
of the valence distributions at larger values of Q , and the value of K. 
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Note that, if higher- twist effects and other uncertainties were negligible, the following 
dependence of K on Q should hold 



K(Q' ) = K(Q ) + - Eq6 '. In % « K(Q ) + 0.770 In 9k 



With #(0.6 GeV) = 1, ([TJD predicts K(2 GeV) = 1.93, to be compared with 2.42, the 
outcome of the high-Qo fit- Two facts account for almost all the discrepancy. First, the 
X 2 analysis suggests that #(0.6 GeV) should slightly exceed unity (K = 1.17) increasing 
the prediction fll8l) to K (2 GeV) = 2.10. Second, the high-Q fit includes only data 
above Q 2 = 4 GeV 2 , while the low-Qo fit includes data down to 0.71 GeV 2 . Indeed, good 
agreement is found if the latter fit is also restricted to the high-Q 2 data. This indicates 
that higher-twist contributions significantly affect at low Q 2 , but also that the leading- 
twist evolution of PDFs is still valid down to Qo = 0.6 GeV. 

Distributions involving the inclusive PDFs of the real photon are not sufficient to disen- 
tangle the non-perturbative part. The high-Qo and low-Qo SaS sets of PDFs describe the 
F2 data equally well [f|], although the non-perturbative parts have very different shapes 
and normalizations. For example, the momentum fractions carried by the perturbative 
and non-perturbative parts are quite different for the two cases but not visible with a real 
photon (P 2 = 0) target: 

Qo = 2.0 GeV : 1 - Z 3 = a em jl.33 + ^ ^ej) In Q " 

Qo = 0.6 GeV : 1 - Z 3 = a em { 0.55 + — ( ]T e 2 ] In 




tt IV 7 0.36 GeV 



in 

2 



U>6 + ^ ? 4 K^} ' (19) 



Additional information can, and has to, be obtained from two sides: perturbative and 
non-perturbative parts lead to differences in the hadronic final states of photon-induced 
reactions ||, and also show a different dependence on the photon virtuality as is evident, 
e.g. from QHp. Further studies, both theoretically and experimentally, to exploit the 
sensitivity to the non-perturbative part are highly desirable. 
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